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Abstract: We consider chiral perturbation theory in a finite volume and in a mixed regime 
of quark masses. We take N[ light quarks near the chiral limit, in the so-called e-regime, 
while the remaining quarks are heavier and in the standard p-regime. We compute in 
this new mixed regime the finite-size scaling of the light meson correlators in the scalar, 
pseudoscalar, vector and axial vector channels. Using the replica method, we easily extend 
our results to the partially quenched theory. With the help of our results, lattice QCD 
simulations with 2+1 flavors can safely investigate pion physics with very light up and 
down quark masses even in the region where the pion's correlation length overcomes the 
size of the space-time lattice. 
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1. Introduction 

It is now becoming feasible to simulate Quantum Chromodynamics (QCD) near the chiral 
limit in lattice gauge theory. Deeper theoretical understanding of chiral symmetry on the 
lattice, development of very efficient algorithms, and a constant progress in the compu- 
tational resources have allowed a reduction of the dynamical quark masses almost to the 
physical point of the u and d quarks. As the quark masses approach zero, however, one 
has to be increasingly careful about finite-volume effects since the correlation lengths of 
the pions, the pseudo-Nambu-Goldstone bosons, diverge in that limit. 

Such an infrared effect due to finite volume can be systematically treated within the 
framework of the low-energy effective theory. Due to the mass gap between the lightest 
particles, the pseudo-Nambu-Goldstone bosons that are generically referred to as "pions", 
and the other hadrons, the heavier particles have an entirely different sensitivity to the 
finite volume. The euclidean partition function of QCD receives contributions from the 
full spectrum, but if the chiral limit is taken at finite volume these higher states give 
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exponentially suppressed contributions. In this way, by varying the volume, one can tune 
to as high accuracy as one wants by including only those degrees of freedom that are 
associated with the pseudo-Nambu-Goldstone bosons. In QCD, one is interested in a 
situation where two of the quarks (the u and the d) are extremely light, while a third (the 
s) is closer to the QCD scale Aqcd> but still light on the scale of the ultraviolet cut-off 
4ttF (where F is the pion decay constant) in the effective low-energy theory. 

It is therefore important to investigate how the finite volume can affect low-energy 
dynamics within the pion effective theory, chiral perturbation theory (ChPT) An 
extreme finite- volume situation is reached in the so-called e-regime where the pion correla- 
tion length l/m n exceeds the size L, T of the 4-dimensional space-time volume V = TL 5 , 

< L,T < — . (1.1) 



Aqcd m* 

The lower bound is to ensure validity of the effective chiral theory, the coupling constants 
of which are the same as those at infinite volume. 

A systematic expansion exists in the e-regime, where all zero-momentum modes of 
the pseudo-Nambu-Goldstone bosons have to be treated exactly. An appropriate power- 
counting in this regime is 

„2 i i T 2 i /rp2 ,2 



- p z ~ 1/L Z ~ 1/T Z ~ e 2 , (1.2) 

in units of the cut-off of the theory. With this counting, the operators in the Chiral 
Lagrangian have different weights than in the ordinary ChPT at infinite volume, known 
as the p-expansion. There is therefore a re-ordering of the perturbative expansion: in 
many cases, the infinite-volume chiral condensate £ and the pion decay constant F (both 
in the massless limit), play a more prominent role than in the conventional large volume 
regime, since next-to-leading order corrections are calculable in terms of these leading- 
order couplings alone. This opens up the possibility of extracting some of these low-energy 
constants from lattice QCD in new ways 

With non-degenerate quark masses one can define an e-regime and a p-regime for 
each of the now mass-split pseudo-Nambu-Goldstone bosons. In particular, one can also 
consider a mixed regime in which some pseudo-Nambu-Goldstones obey the condition for 
the e-regime, while others fall into the p-regime H. For the latter, the counting rules are 
the usual ones of chiral perturbation theory, 

m q ~ml ~ p 2 ~ 1/L 2 - 1/T 2 ~ e 2 . (1.3) 

A typical situation could be u and d quarks so light that the physical pions are in the 
e-regime, but the strange quark mass is such that the physical kaons are in the p-regime. 
Further possibilities open up when one considers partially quenched theories. In those cases 
one can imagine situations in which all physical u, d and s quarks are in the p-regime, while 
valence quarks corresponding to all or some of these are taken closer towards the chiral 
limit, and thus end up in the e-regime. Such situations could perhaps be realized in the 
context of mixed-action lattice simulations where dynamical configurations are generated 
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with physical quarks that are in the p-regime and can be well treated by, say, ordinary 
Wilson fermion actions. Valence quarks, which are taken to the chiral limit, could then be 
of, say, overlap type. Another situation could be the use of overlap quarks that are all or 
partly in the p-regime, while also overlap valence quarks are taken to the e-regime. 

It is our belief that all these possibilities must be and will be explored in future lattice 
gauge theory studies. The present paper gives an analytical formalism for studying cor- 
relation functions of pseudo-Nambu-Goldstone bosons in that setting. Apart from quark 
masses, the two remaining limitations are the approach to the continuum limit, and the 
finite volume. In the present framework the finite volume V S> 1/Aq CD is used as a tun- 
able parameter with which to extract physical observables. The only extrapolation needed 
will thus be the one associated with taking the continuum scaling limit. With the new 
analytical formulas for partially quenched correlation functions available one can extract 
far more information for a given number of lattice configurations. In addition, with the 
mixed-regime predictions we also provide here one is effectively covering the full SU (3) 
flavor sector of QCD at low energy. The extent to which the s quark at the physical point 
is light enough to provide a good description in terms of chiral dynamics remains to be 
tested in detail on the lattice. 

One important feature of the e-regime is the strong sensitivity to the topology of the 
gauge fields, a direct consequence of the finite volume. A crucial ingredient in making 
the e-regime so useful for lattice gauge theory computations is the fact that the zero- 
momentum integrals can be performed analytically at fixed topology, typically in terms of 
Bessel functions. With different detailed predictions for each sector of topological charge 
there is a wealth of analytical results that can be used to confront numerical lattice data. 

These analytical predictions are non-perturbative in the gauge theory coupling. Re- 
markably, some of the leading-order e-regime results were first derived on the basis of chiral 
Random Matrix Theory j7||. It is particularly simple to derive analytical expressions for 
Dirac operator eigenvalues in that formulation ||, an d it is well understood how to go 
between the two formulations ||. These features all carry over into the mixed regime. 

The computation of meson correlators in the e-regime was first performed in m and 



extended to both quenched and unquenched QCD at fixed topology a few years ago [10] - 
fj~2[l . The effect of a coupling to isospin chemical potential has also been considered [13]. 
Partially quenched ChPT (PQChPT) in the e-regime was done for the chiral condensate 
itself in [14]. Recently, partially quenched space-time correlation functions were computed 



in several channels (pseudoscalar, scalar, and the left-current) of meson correlators ]15|, |(|. 
Also, the first computation of three-point correlation functions relevant for weak decays 



in the e-regime was done in [16]. Some of these studies have led to determinations of 
the leading order low-energy coefficients S and F 1 at various sectors of fixed topology in 
quenched lattice simulations (see, e.g., 18]-p8|), as well as the low-energy couplings of 



the AS = 1 Hamiltonian from three-point functions p9| , 30]. Difficulties associated with 
the quenched approximation have been discussed in |l(], [l^] . Recently, several groups have 
successfully extended this to full QCD in or close to the e-regime. This has been done both 



1 See, e.g. ref. [jfjj for a recent summary of results. 
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on the basis of Dirac operator eigenvalues and space-time correlation functions [31]- [40]. 

In this paper, we present the results for various meson correlators in the mixed regime 
of ChPT where iVj light quarks are in the e-regime while = Nj — Ni quarks remain 
relatively heavy and belong to the standard p-regime. We have used two different methods 
to treat this regime. The first uses the same mixed-regime perturbative expansion that was 
introduced in Q. As a check, we have also used a new perturbative approach which has the 
advantage that it provides a smooth interpolation between the e-regime and the p-regime. 
The expected matching between the mixed regime and the standard e-regime is trivial to 
check in that formalism. The two methods should agree to all orders, and we have checked 
explicitly that they do agree at least up next-to-leading order in the mixed-regime power 
counting. 

We treat the most general non-degenerate case where the required non-perturbative 
zero-mode integrals are performed according to Ref. [15]. The two-point functions of the 
light sector, for the pseudoscalar, scalar, axial, vector channels, are then computed. They 
can be used to extract the leading low-energy constants £ and F. Because we work at 
next-to-leading order, there is also explicit dependence on some of the L^s. In principle 
these low-energy constants can be determined from fits to varying quark masses in the 
heavier sector, as will become clear below. The pseudoscalar and scalar channels for the 
disconnected diagrams are also given. We easily extend our results to the partially (and 
fully) quenched theory by applying the replica method. One can confirm that our formulae 
reduce to all previously derived limiting cases of both the degenerate iVj-flavor theories 
and the fully quenched theory. There are new isospin-breaking effects when the u and d 
quark masses are split, and the existence of these terms can be used to extract additional 
information from the correlators. The new more general expressions should be helpful for 
future lattice gauge theory simulations that aim at approaching the chiral limit. 

We start in Section ^ by reviewing the mixed-regime perturbative expansion of || . The 
results for the two-point functions at next-to-leading order are presented in Section ||[ An 
alternative new approach is also briefly described there. The calculations are completed 
in Section || by explicitly performing the zero-mode integrals for the full, the partially 
quenched, and the fully quenched theories. As a check on our results, we note the complete 
agreement between our two approaches and the correct matching between the mixed and 
pure e regimes is then also explicitly confirmed. In Section [5], we give an explicit example 
for Nf = 2 + 1 theory presenting the pseudoscalar and axial vector correlators. Conclusions 
and an outlook for the future are presented in Section 0. 



2. Chiral Perturbation Theory in the mixed-regime 

In this section we review the perturbative expansion of the chiral Lagrangian that was 
introduced in || to treat the mixed regime. It incorporates features of both e and p 
expansions, allowing for the simultaneous presence of quarks with masses corresponding 
to these two regimes. The existence of such a mixed expansion will be useful for lattice 
simulations at the physical points of the three lightest quark flavors u, d and s, or, more 
modestly, simulations where only associated valence quarks are taken to that limit. 
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Let us consider an iVy-flavor theory in a finite volume V = L 3 T, 

F 2 £ 
C = —T,[d^U{x)^d^U(x)] - -Tr[M^U(x)U e + ulu(x?M\ +••-, (2.1) 

where U(x) £ SU(Nf) and £/# = exp(i0/Nf)l. Here is a QCD vacuum angle, introduced 
here only in order to be able to project on fixed gauge field topology by doing a Fourier 
transform in 9. As usual, £ denotes the infinite- volume chiral condensate in the massless 
limit, and F is similarly the pion decay constant in the chiral limit. Note that there 
are next-to-leading order terms, indicated here by ellipses, each of which correspond to 
additional low-energy constants denoted by, in the SU(3) case, L,. 

For the mass matrix M = diag(mi, ni2 • • •), we consider the most general non-degenerate 
case, where we have Ni light quark masses in the e-regime: 

M hh =m h ~0(l/V) 1 (2.2) 

while the other = Nf — Ni quarks are heavier: 

Mhiht = m hi ~ Oil/V 1 / 2 ), (2.3) 

in units of the cut-off of the theory. Here and in the following, we put a subscript I for the 
light sector and h for the heavier sector and denote the mass matrices in those sectors by: 

M t = PiMPi M h = P h MP h , (2.4) 

where Pi,Ph are projectors on the light and heavier sectors respectively. The working 
assumption is of course always that chiral perturbation theory is meaningful even for the 
heavier sector. 

In ref. || an expansion was proposed according to the following counting rules: 

ZV~0(e), L,T~0(l/e), Mu ~ C(e 4 ), M hh ~ 0(e 2 ). (2.5) 

An inspection of the pion propagator shows that the zero modes of the Nambu-Goldstone 
fields associated with the generators in SU(Ni) need to be treated non-perturbatively. 
All the remaining zero-modes are perturbative when Ni,Nh ^ 0. Some subtleties appear 
however in the partially-quenched case where all the light quarks are quenched, that is the 
replica limit Ni — » 2 . In this case, it is easy to see that the Goldstone field associated to 
the generator T v , 

T,= (^f l ?A (2.6) 

gets massless in the replica limit Ni = 0. Here 1;/^ are the identity matrices in the light 
and heavier sectors. Note that T v looks ill-defined when Ni = but keeping Ni finite until 
the very end of the calculation, one sees that the replica limit Ni — > can be safely taken. 



2 As usual the fully-quenched case Ni + Nu — requires the presence of the singlet to be well-defined, 
but as long as iV; + Nh 7^ the singlet decouples. 
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To treat all cases on the same footing, we therefore consider the following parametriza- 
tion 3 : 

t /(x)=exp(^)(^ I °Jexp( j ,T,), (2.7) 

where Uq £ SU(Ni) is a constant matrix, r\ is the zero-mode of the Nambu-Goldstone field 
associated with the T v generator. The ^ fields contain the non-zero modes corresponding 
to all Nambu-Goldstone fields, and also all zero modes of those degrees of freedom that are 
not treated separately. They therefore satisfy the constraints 



J d A x Ti[T a i{x)} = J d^x Tr[T v £(x)] = 0, (2.f 



where T a is a generator of the subgroup SU{N{). Note that the zero-mode of the T v 
generator is not included in the £ field (it is projected out by the second constraint in 



eq. ( p.8[) ), and included explicitely in the last term of eq. (2.7). 

We are interested in computing the correlation functions in sectors of fixed topology. 
Following the same derivation in § we rewrite 

ie 

i e 2N i T; D i 

U(x)U e = U(x) 1 J_ = exp 



>2Af h I 



h 




where we have defined 



ri — - ri + 9 , 

and Uq 6 U(N[) with det Uq = e^det(CTo) = e l6 '. The partition functional in sectors of 
fixed topology can then be written as: 

Z v ~ J [d£\[dfl\j uN [dU ] J(0 det(U ) u exp (- J d 4 x£(£,fi,U ) 



(2.11) 



where as in the standard e-regime, the projection on fixed topology results in the enlarge- 
ment of the zero-mode integration from SU(Ni) to U(Ni). </(£) is the Jacobian of the 
change of variables of eq. ( |2.7D , According to the power-counting of eq. ( |2.5| ), it can be 
shown that a consistent power-counting for the fields ^ is: 

e~0(e), (2.12) 

therefore both the Lagrangian and the Jacobian can be perturbatively expanded in powers 
of £. At next-to- leading order we find [|, |: 

J(0 = 1 - ^ / d'x Y, Tr ^ 2 " ( T ^) 2 K^) + 0(e% (2.13) 

atSUiN^UTv 



3 In reference |H| a different parametrization was considered for the case when some or all of the light 
quarks are dynamical. In that case the rj zero-mode is also perturbative and can be included in £. It turns 



out that the parametrization of eq. (2.7) simplifies the calculations and allows one to consider the full and 



partially quenched cases on the same footing. Therefore we consider only eq. (2.7) in the present paper. 
We have checked that both give the same result in the full case. 
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The Lagrangian can also be obtained as an expansion in e: 

£ = £ (4) +£ (6) + > 



(2.14) 



with terms up to 0(e 4 ), up to 0(e e ), etc . Concerning the integration over the variable 77, 



we can perform a saddle-point approximation following the derivation of 
order Lagrangian is found to be: 



]. The leading- 



= Tr [d^] - -Tr \Mi(U + 0}) 



-I Tr 



M h [i 



Fjj_. 
~2Nh' 



v _ 



(2.15) 



This quadratic form implies also a power-counting of f) ~ e. According to this rule, the last 
term in eq. ( 2.1 5| ) could be treated as a perturbation. This is true as long as v ~ 0(e°), 
as is usually the case in the e-regime. However, in the partially-quenched case iV; = 0, 
the distribution of topological charge is controlled by the heavy quarks only. Indeed the v 
dependence of the leading-order partition function is found to be 

/ ,,2 



Z, 



LO 



oc exp 



1 



'U(Ni) 



[dU ] det(U y exp -Tr Mi(U + U 



which in the case A^; 



implies: 

W 2 ) 



-VF 



(2.16) 



(2.17) 



V -4, 

a scaling that makes the last term in eq. ( 2.15| ) of C(e 4 ), and therefore of leading-order. In 
order to recover the results at 9 = by averaging over topology, it is therefore necessary 
to keep the last term in eq. ( |2.15| ) in the leading-order Lagrangian, or equivalently assume 
that v ~ e _1 . This is not necessary however as long as N\ > 0, since the distribution of 
topological charge in that case is controlled by the light quarks. 

It is straightforward to derive the propagator for the £ fields in the light or mixed 
sectors from eq. ( |2.15 ), validating the power-counting of eq. ( [2.12 ) and ensuring that the 
replica limit iVj = is well-defined: 

8 hh 8i 2 i 3 A(x - y, 0) - 5 h i 2 5 hh G{x - y, 0, 0) 



filial) 6s J* (l/) 

finite) €h 2 h(y) 
fiii 2 te)&iift 2 (y) 



2 L 
1 



dhhShih 2 A 



M 2 hlhl 



-5 hh 5 hlh2 G (x - y,0, Ml lhi ) , 



(2.18) 
(2.19) 
(2.20) 



while in the heavy sector there always appear the combination: 



611/12 te) 



Ffj 
2N~h 



Ffj 
2N~h 



8h... 



S hlh4 S h2ha A(x - y,Ml lh2 ) - 5 hlh2 S h3h4 (G(x - y,M% lhl ,M% 3h3 ) 
+G (Ml hv Ml h3 )) 



(2.21) 
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where 



1 ptpx i 

^ M2 )^V^WTW^ A(z,M 2 )^A(*,M 2 ) + — (2.22) 
G(x, Mf, Mi) = -Y —, r (2.23) 



t^o (P 2 + Ml)(pi + M*) ^ + E 



'h p 2 +M 



2v l 



and 



Ml hl = {m h + m hl )^, (2.25) 

is the mass of the meson fields made of the heavier quarks. The summation Ylp^o * s taken 
over the 4-momentum 

p = 2Tr(n t /T,n x /L,n y /L,n z /L), (2.26) 

with integers n^'s. Note that the term Go is formally of higher order if v ~ 0(1). 

In this work, one encounters G{x, M 2 , M 2 ) with Mi = M2 = only, both in full and 
partially quenched theory 4 . In the full theory, G(x,0, 0) can in principle be rewritten in 
terms of A's, as one would have expected on general grounds. In the case of N[ = 2 and 
Nh = 1, which is the phenomenologically most interesting case, for example, 

1 pipx 1 1_2 

G(x, 0,0) = V J2 a(2 ■ I ^ = 5 A (*> °) " 6 A(X ' 3 M ' ft) - (2 - 27) 



But we keep using the notation of G(x, 0, 0) for simplicity in both the general case with 
Ni + Nh flavors and the partially quenched case, where a double pole appears. 

Correlation functions are obtained by inserting appropriate source operators in the 
above partition function and taking suitable functional derivatives ||. The U(Ni) integral 
over zero modes Uq is then done exactly, while the £, fj integrals are treated perturbatively. 
We return to the zero-mode integrations in Section ^. Here we will be working at next- 
to- leading order in the perturbative expansion and therefore up to the term in the 
Lagrangian contributes: 

£(6) = [(^(*)^)) 2 - («*)) 2 e 2 o«o] - ^i Tr [-^ 4 (*)] 

y r / - Nil fiV^ T 

+ -^Tr [Mi (e(x)Uo + Ufe 2 (x))\ + — ^Tr [M h ] Tr [d^d^x)] 

- 16^Tr [M h ] Tr \Mi (u + u£)} + . . . (2.28) 



F 4 

where the ellipses indicate terms of the same order that involve only £^ or an d do not 
contribute to the observables where valence quarks are only in the light sector, as we will 
be considering in this paper. 



4 The fully quenched case needs special care; it will be discussed later. 
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3. Two-point correlation functions 

Correlation functions are obtained by inserting appropriate source operators in the above 
partition function and taking suitable functional derivatives Q. Since we consider the 
fully non-degenerate theory, we have to treat all possible Nf x Nf bilinear quark operators 
separately. We therefore define 

P v {x) = iqi(x)-/ 5 qj(x), S tJ (x) = qi(x)qj(x), 
A l i{x) =iq i (x)^ li qj(x), Vff \x) = iqi(x)jpqj(x). (3.1) 

The corresponding operators in ChPT are to the leading order given by 

P ij (x) = (\U{x)U e ] i3 - > ( 3 - 2 ) 

S ij (x) = | {[U{x)U e ] tJ + PlU^x)}^ , (3.3) 

= i^-[d»U(x)U\x) - d^(x)U(x)} l3 , (3.4) 

V?(x) = i^[d^U{x)U\x) + d^(x)U(x)] ir (3.5) 

The conventional irreducible representations are obtained by appropriate combinations of 
i's and j's. The charged pion-type meson operator and the neutral one are, for example, 
given by (we simply denote 1 for the up quark and 2 for the down quark) 

P^(x) = ^(P 12 (x) + P 21 (x)), and P*\x) = ^(P n (x) - P 22 {x)). 

(3.6) 

In the following, we use indices v and v' in order to specify the valence sector which in this 
paper is always taken to be in the e-regime. 

In Figures [l] and |2| we show the Feynman diagrams resulting from the £ integration 
that contribute to the current and scalar propagators at next-to-leading order in the e- 
expansion. The scalar correlators start at 0(e°), while the first contribution to the currents 
is C(e 2 ). Note that disconnected diagrams contribute because they are connected through 
the zero-mode integrations. We also assume here that the operators are separated from 
each other and the usual contact terms are not included. 

For the practical purpose of comparing to lattice QCD simulations, we will present re- 
sults in the non-singlet irreducible representation, or " charged-pion" type correlation func- 
tions, with the zero-momentum projection (integration over 3-dimensional space), namely, 



KAt) = \J d 3 x((P^(x)+P^(x))(P^'(0)+P^(0))) u ^ (3.7) 

S c vv ,(t) = lj d 3 x((S™'(x) + S^(x))(S™'(0) + ^(0))}^, (3.8) 

A%A*) = \J d 3 x((Af(x)+A^(x))(Af(0)+A^(0))) u ^ (3.9) 

KAt) = \J d 3 x((V^'(x) + V v ' v (x))(V o m '(0) + ^(0)))^, (3.10) 
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Figure 1: Diagrams contributing to the scalar and pseudoscalar connected correlators. The lines 
are £ propagators. Squares indicate the scalar and pseudoscalar operators. The filled dots indicate 
a mass insertion from the Lagrangian. Empty dots indicate the insertion of an operator coming 
from the NLO Lagrangian, and they are also labeled with the subindex of the associated coupling 
constant, L t . 





Figure 2: Diagrams contributing to the vector and axial correlators. The lines are £ propagators. 
Squares indicate the vector or axial vector operators. The filled dots indicate a mass insertion. A 
cross indicate a contribution from the Jacobian. Empty dots indicate the insertion of an operator 
coming from the NLO Lagrangian, and they are also labeled with the subindex of the associated 
coupling constant, Lj. 

where we assume w/i/. 

We also present the "disconnected" contributions for the scalar and pseudoscalar s, 

VUt) = j d 3 x{P™(x)P v ' v '(0))u£, (3.11) 
SUt) ee J d 3 x(S vv (x)S^'(0))u^ (3.12) 

which are useful to estimate the finite size contributions from the chiral fields to the r( 
meson correlators. 

To simplify the t-dependence of our expressions, let us define 

M*/r) = ^J d 3 xA(x,0) = \ (| - £) 2 - 1, (3.13) 



and 



•(*) ee J d 3 x G(x,0,0). (3.14) 
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Note that the latter still depends on the heavier quark masses. In the appendix |A|, we 
list how to perform the zero-momentum projection of various combination of A(x, M 2 )'s. 
Defining also 



S _ 1 
S = 1 ~ F 2 



F _ 1 
F = 1 ~ 2F 2 



iV, A(0, 0) + ^ A(0, M 2 hh /2) - (5(0, 0, 0) - 16L 6 £ Mj 
ft 

JVi A(0, 0) + £ A(0, M^/2) - 8L 4 £ M ( 



(3.15) 
(3.16) 



and 



(3.17) 



the results for the pseudoscalar and scalar non-singlet (connected) correlators can be writ- 
ten as 



-L 6 



K 



0(Ni) 



+ 



S 2 



S c vv ,(t)= L 3 ^- 



K 



0(Ni) 



NLO 2F 2 

T 2 



NLO 2F 2 



K}± Nl) Th x ( - ) +K?} in >r{t) 



-O(iVj) 



(3.18) 
(3.19) 



and their disconnected correlators are given by 

2T 2 



Vi,{t) = -L 3 t 2 
Si,(t)= L 3 t 2 



K 



2W) 



+ 



K 



mi) 



NLO F 2 

2Y 2 



NLO F 2 



K*W)Th x 



1) 


_ K -2iN l) 




1) 









(3.20) 
(3.21) 



The tC functions represent the zero- mode integrals over U{N{) and they depend only on 
the light quark masses (to simplify the notation we denote Uq by U in this section): 



kT i) (M) = ^((U m> + U v , v ± ul, v ± ul,f) um 



1 



rt \2\ 



(3.22) 



ZCi (7Vi) ({ W }) = 1 ± \{U vv U v , vl + Ulul, v ,) um ± i(C7^, + £/ 2 ,„ + ^.c.)^), (3.23) 



/c! (7Vi) ({M) = ± t4)(£W ± ul >v ,)} um , 



where averages are over zero modes: 

((■■■))u m - f dU(...)(detUye^^ u + u ^. 

JU(Ni) 



(3.24) 
(3.25) 

(3.26) 



The label Ojvlo i m pli es that the integral must be computed with £ instead of S. We will 
present the explicit results for these integrals in Section |j. 
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For the axial and vector (connected) current correlators we obtain: 



?2 r 



A c vv ,(t) 



2T 
S 



j: 



0^) 



T 

NLO 2V 



N h 



i( Nl ) 2Nj_ / ( Nl ) _ ( Nl ) 



iV^ 00 + ^^ (M h \/2) 
Thi 



(3.27) 



?2 r 



v c vv ,(t) 



2T 



T 

NLO 2V 



N h 



A^ 00 + J>^ (M 2 ,/2) 

h 



j: 



0(N t ) 



Thi - 



where we have defined 



^ (M 2 )etAa(0,M 2 ) 



T 2 d - 

— &oo = T— A(0,0), 



and the functions are given by: 

s^HM) = l(Uw + ul) um , 
J± Nl) (M) = i ± ( Uvv ' u v v ' + + h - c -^w) 

4 W) ({ W }) = ((2rrv±m v ),S^ 

{U$ v (UMiU) vlv > + U^UMiU)^ +h.c.) u{Nf) 



(3.28) 
(3.29) 

(3.30) 
(3.31) 

± (u <-> t/) (3.32) 



We stress that all the heavier mass dependence is explicit in the results of eqs. ( |3.18 )- 
(1319D , (g!20|) - (^2lD and (|3i27|) - (^28D since the zero-mode integrals involve the light sector 
only. We also note that these results agree with those obtained for the special case of the 
left-handed current two-point function obtained in |]]. 

Next, we need to discuss the ultraviolet divergences of A(0,M 2 )'s and similar ones 
associated with G"s. The explicit form in finite volume is given by Q, 

M 2 



A(0,M 2 



16vr 2 



(lnM 2 + ci)+ 5 i(M 2 ), 



(3.33) 



where c\ represents the logarithmic divergence which is independent of M and the volume, 
and g\ denotes the finite volume correction . The numerical evaluation of g\ is discussed 
in Section || 

Since F and m^S are not renormalized at infinite volume, the logarithmic divergence 
c\ must be absorbed in a renormalization of L^s. The mass-independent shift in 



L4 — > L4 + 



(16vr) : 



r Ci, Lq — » Lq + 



(16vr) 5 



* + 4 



Cl 



(3.34) 



is enough to give finite results in the above correlators. This shift is exactly the same as 
at infinite volume 111. 
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Note that the formally divergent expressions other than A(0,M 2 ) (for non-zero M) 

/2' 

-1 



A(0,0) = --^, (3.35) 



"°° = ,^ w 4smh4|T/2) + ^' (3 - 36) 

become finite after dimensional regularization. 

Finally we note that the dependence on the heavier quark masses is as expected on 
general grounds (see also the discussion in Q). Indeed, up to exponentially suppressed 
finite- volume corrections in M^L, the correlators above coincide with those in the e-regime 
for Ni light quarks as if there were no heavier quarks whatsoever. The only remnant of 
the heavier quarks is seen in the modified low-energy couplings £ and F, i.e. by the terms 
that depend on M^h in £ and F. This is as usual in chiral perturbation theory. 

3.1 An alternative mixed-regime expansion 

As a check on our results, we have performed the same calculation by means of an al- 
ternative method where the parametrization of fields is as in the standard e-regime. The 
counting rule we use, however, is the same as the one in the standard p-regime for the heavy 
flavors. All zero modes in the full SU(Nf) group are then treated non-perturbatively. Such 
a parametrization has the advantage that the matching to the e regime is smooth by con- 
struction. 

The result of this alternative scheme leads to definitions of £ and F which are identical 
to eqs. fl3.15|) and (|3.16| ) except for the replacements A —* A and G —* G. Similarly, all 
other results presented above are reproduced with the only difference that now all zero- 
mode integrals are performed over the whole U(Nf) group and therefore depend on all the 
quark masses, including the heavier ones. 

In contrast to the results presented in Eqs. . (gl8|) - (|3.21j) and (^27|) - (|3.28D ), in this 



alternative approach one can take the limit M^h — * smoothly. The results then coincide 
with those fully in the e-regime. Indeed, our results in that limit agree with partially 
quenched scalar and pseudoscalar correlators for non-degenerate masses that can be found 
in ref. [15]. The left-handed current correlator can be found in f|, and our present results 



also reproduce that special case. 

The reason that the matching limit is smooth in this parametrization is because the 
zero-momentum modes of the massive mesons are resummed, while in the expansion of 
Section [2], they are treated perturbatively. The two results should therefore coincide when 
the zero-mode integrals of the [/(iVy)-theory are expanded to the appropriate order in 
l/Hh ~ C(e 2 )- I n the next section we will show that this is indeed the case. This provides 
a rather non-trivial consistency check on our results, and it confirms that the expected 
matching between e and mixed regimes actually holds. 
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4. Non-perturbative zero-mode integrals 



In this section, we complete the calculations of the correlators by giving explicitly the zero- 
mode integrals defined in eqs. ( p . 22| )- ( p725 ) and ( 3.30| )- ( |3.32| ), in the full (unquenched) , 



the partially quenched, and the fully quenched theories. Here we present the results of 
general partially quenched calculations. As is well-known, the results for the full theory 
can be viewed as special cases, obtained by equating the valence quark masses to those of 



the sea quarks. The essential ingredient is the functional [42], 



7 v ( r v, det[/i| 1 Ju+j-l(Vi)]i,j=l,-n+rn 

nMM) ^^(^ _ ^ U n+ r=n+i(f ,2 _ ft > ^ 

where /ij = m^V. Here ^7's are defined as J' u+ j-i(fj,i) = (—iy~ 1 K u+ j-i(fii) for i = 1, ■ ■ ■ n 
and Jv+j-i{[ii) = Iu+j-iim) for i = n + 1, • • • n + m, where I u and K v are the modified 
Bessel functions, m = N v + N[ denotes the number of (light) quarks, from which N v valence 
quarks are quenched by the n = N v bosonic quark contents. Since we are interested in 
mesonic two-point functions, we need n = N v = 2 at most. 

From this functional, by taking appropriate derivatives with respect to the parameters 
Hi, one can derive all the required integrals, both in the full as in the partially-quenched 
limits. The technical steps of our calculation have followed those of Ref.jL^] and in this 
section we simply show the final results. Details of how this can be used to compute all 
relevant group integrals are presented in Appendix [B|. An important relation can be derived 
from Ward-Takahashi identities (see Appendix |C]) that holds for the full, partially-quenched 
and quenched cases: 

jl™ = 2K ± ™ V/ )(SW ± S™) T - (4.2) 
As building blocks, let us define two quantities, 
x£ Q (^„,{M) _ d 

= - hm — lnZ 11+N {ii b \fj, v ,{fi s }), (4.3) 

hm ' ' — — — — . (4.4) 

where the sea quark mass dependence is denoted by {[J>s} = {ft, m N[ }. 

Let us also give their analogous expressions in the unquenched theory (we need only the 
case where the valence mass is equal to one of the light sector, m v = mi in the e-regime.), 

v = ~— lnZ Ni (W = lim v > ( 4 - 5 ) 



nfull/ r "n — ^Wi^W2-^JVj ({/*«}) ,. r>PQ/ r T\ 



(4.6) 
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and the fully quenched limits, 



lim E^ Q (w{w}) = T,l Q (ii v ) ^ ^ 

v 2 

lim £)P Q (/i vl ,^2,K}) = 1 + • (4.8) 

{Ms}— >oo WlW2 

4.1 Explicit results 

With the above expressions, one can calculate all the non-perturbative integrals we need 
to evaluate J± etc. Further details can be found in Appendix |b|. 

1. Full (unquenched) theory 

We start by listing the results for the full (unquenched) theory, where the valence 
masses are equal to those of the sea quarks (m v = mi and m v i = my). 

5," V ''fa,fa,}) = S " ( ^ {/ " }) , (4.9) 

K ™>(„, w „ W ) = -J*- ( iPfa.W) T ^WT) , (4 , 10) 

Mi T W V ^ ^ / 

/Ci W) ( W , W ,{/i s }) = 1± ( D ™\ f , l , fM ,,{^}) + —) , (4.11) 

4 ( ^Vz,W',{M) = ^ ull (w,W,{w}), (4-12) 

/C 2 _ W) ( W ,W',{w}) = — , (4.13) 

WW 

w - w V s s / 



J-° W) ( W ,W,{W» = 1± (^(MJ.Mi'.W)-— 1 , (4-15) 

V ww / 

J ± v ,; (/4,W>{M) = 2(mi ±m F ) I ± 



T f£ - J° m ) , (4.16) 
where = {fx h , m 2 , ■ ■ • /^}. 
2. Partially quenched theory (iVj 7^ 0) 

The partially quenched results where the valence masses are different from the seq 
quark masses, are obtained analogously for the case JVj 7^ 0, 

4 Nl \»vM) = EPv \ M \ (4.17) 

= / ^ Q (w{wl) T ^ Q (w,{w}) \ (418) 
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K}^ Nl \^ v , /v,{ M) = 1± ( D^ Q (/^,/v,{Ms}) + 



£ + ( "Otx«,/v,{M«}) = ^ Q (^,^',{/"s})> 
K. 2 J: Nl \fj, v ,iJL V f,{fj, s }) 

/C«(/x„,/v,{M) 



1 



(4.19) 

(4.20) 
(4.21) 



2 ,,2~ I Mv' ^ M» ^ |l) 4 -^J 



MS - Mi' 



^ W) (^,/v,{M) = i± (^ q (m,,/v,{M) 



J7"2 (m*j,/V,{Ms}) = 2(m„±m„/ 



S^ Q (/x„,{/i s }) £^ Q (/v,{M) 



3. Partially quenched theory (iV; = 0) 



(4.23) 



(4.24) 



In the case with JVj = 0, one needs the /uZZy quenched integral over Z7o: 

^ Q (/i.) 



i-0(0)/ \ 

yv ± (fj, v ,fi v >) 



/c 3 (°)(m„,^) 



±2 



A*i> ~F Mt 



/Cj_ (0) (/^;/v) = 1 + 



H"t)H"u' 



^ Md ^ 



u 2_„2 M,' 

Md /v y 

0(0), 



J± (Vv, Vv>) = 2(m v ±m v >) 



Ey Q (/i„) ££ w (/v) 



FQ, 



(4.25) 
(4.26) 

(4.27) 

(4.28) 

(4.29) 
(4.30) 
(4.31) 



4. Fully quenched theory (iV; = = 0) 



When iVj = Nh = or the theory is fully quenched, the zero-mode integrals we use 
are the same as the partially quenched case with jtyj = above. But we need further 
to include the singlet degree of freedom for the non-zero modes, with additional low- 



energy constants a and tuq as quenched artifacts [43, 44, 11, 12 1. This results in the 
modification of r(t) to 

,2r 



24 



T J \f~ ) ~ 30 



+ aT/n(t) , 



(4.32) 
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where N c denotes the number of colors. See e.g. ref. 11] for details. 



Using the unitarity condition given in ref. [15], we have checked that all of the above 
expressions precisely reproduce the known results for degenerate iVj 7^ flavors (setting 
Nh = 0), and the quenched results (the iVj = Nh = limit), obtained earlier jLl], Q. 

4.2 Equivalence of the two mixed-regime expansions 

The apparent difference between results obtained in the two mixed-regime expansions con- 
sidered in section ^ arises from the contribution of the zero-momentum modes of the heavy 
mesons. They are computed perturbatively in the first case and resummed in the second. 
In order for the two results to agree, an additional expansion in ~ e 2 of the zero- mode 
integrals for U(Nf) must of course be performed so that only terms at subleading order 
in the e expansion are consistently kept in the correlators. Performing this expansion one 
finds 5 : 

j°J Nf) m, w) = Jl m ({*«}) (i-EttVe < 

V h ^ h J h ^ h 

jl w m,{» h }) = jl m ({»i})+o(±- 
4 W m, M) = 4 W) m) f 1 - E - )+ ° (-) 2 - 

Kl (Nf) ({fit}, {fi h }) = JCl {Nl) ({/I,}) + f — ) , n = 1, 2, 3. (4.33) 



.MA 

We have here denoted 



Mi 



ttfi-Efl- (434) 



Using these expansions the results from the two different schemes agree. 

As another non-trivial check in the opposite direction, one can also confirm that a 
fully perturbative approach as in the p-regime, where all of Nf = Ni + iVj, flavors are 
perturbatively treated, is consistent with our results in an unrealistic limit FL S> 1 while 
M n L < 1 kept. 



5. Useful examples for 2+1 flavor theory 

In this section we give some explicit examples that are useful when comparing with lattice 
QCD simulations. Here we consider the 2+1 flavor theory where the up and down quark 
masses are degenerate, m u = and different from the strange quark mass m s . We choose 
the low-energy constants to be the phenomenologically reasonable values F = 90MeV, 
S 1 / 3 = 250MeV, L^(0.77GeV) = 0.1 x 1(T 3 and Lg(0.77GeV) = 0.05 x 10~ 3 . 

5 We have checked these expansions in several special cases with a rather small number of flavors. 
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For the calculation of g\ (M 2 ), we use an expansion in the modified Bessel function 



II, 



M 



\ni\<n m c 

Yl 4vr 2 



■Ki(M\a\), 



(5.1) 



where the summation is taken over 4-dimensional vector = (noT,nxL,n2L,nsL) with 
integers n^'s. Truncation above at n max = 5 already shows a good convergence when M > 
200 MeV and L = T/2 = 2 fm, for example. 

In this theory, and for the cases we will consider, one can express G(x, 0, 0) in terms 
of A(x,M 2 ): 



G(x,0,0) = - [AA(x,M*)+BA(x,0)+Cd M 2A(x,0)} , 

and, therefore, 

r(t) = / d 3 x G(x, 0, 0) 



(5.2) 



A 



1 



cosh(M^(t - T/2)) 
2M V sinh(Af„r/2) ~ M 2 T 



+ BThi(t/T) + CT 3 h 2 (t/T) 



BThx{t/T) + CT 3 h 2 (t/T) 



A 



Mr,t\ 



+ 0(e" M " 



(5.3) 



where 



M*/r) = 24 



TJ \T ) 30' 



(5.4) 



and A, B, C, M„ are functions of the p-regime masses only. The term proportional to C 
only appears in the case with Ni = 0. With this set up, one obtains 



A(0,0) 



01 



<9 M2 A(0,M 2 )| M2=0 



J^ln/x^ 1 / 2 -^, (5.5) 



where (3± and /?2 are the usual shape coefficients and fi su b (=0.77GeV in this section) is the 
subtraction scale. 

With this input, one can now calculate meson correlators on the basis of our expressions. 
In the following, we will give two examples where in both we let the volume size be given 
by L = 2 fm. One is the case where the physical up and down quarks are in the e-regime, 
i.e., Ni = 2 and Nh = 1. The other is the case with rather heavy sea quark masses, i.e., 
Nh = 3, but the valence quarks are taken to the e-regime. 

As seen below, the 1-loop corrections to the condensate and decay constant are con- 
siderable even in the limit V — > oo because of large strange quark mass. Recently, it has 
been argued that Nf =2+1 flavor ChPT at NLO may have difficulty in fitting lattice QCD 
data p6|, |47|| . It is clearly important to check whether NNLO contributions are essential 
for analyzing such lattice results at the physical s-quark mass, or if the strange quark is 
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simply out of the region where ChPT provides a useful expansion. If the latter case is true, 
one would need to integrate the strange quark out and use an "effective" Nf = 2 ChPT. In 
this paper, we do not wish to address this issue and hence just give the NLO formulae for 
the Nf =2+1 theory. Even in the Nf = 2 theory one may be interested in keeping the u 
and d quarks in the p-regime, while taking the valence quarks masses to the e-regime. Our 
formulas given in this paper easily extend to that case, but we do not display them here. 

5.1 The case with N t = 2, N h = 1 

Let us first choose m u = = 2 MeV, m s = 110 MeV, where the physical pions are 
certainly in the e-regime in a volume as small as L = 2 fm. 
In this case, the coefficients of eq. ([O) are 



A 



D 



(7 = 0, r{t) 



1 1 

-ThAt/T) H tt- . 

2 w 1 6M 2 T' 



(5.6) 



where M 2 



3 lvl si 

are then given by 



E 
£ 

F 
F 



1 



1 



1 

1 

2F 2 



= 4 ™pP ■ The 1-loop corrections to the condensate and decay constant 
3A , M 2 K l M 2 K Ml . Ml 1 



+ 



2VV 

y/V ' 16vr 2 



16vr 2 



f^sub 



+ 



M 2 M 2 



^lub 



96vr 2 
16Ll(n sub )Ml 



6M 2 V 



32L r 6 (fi sub )Ml 



(5.7) 



where Ml 



S/F 2 , and we have neglected exponentially small gi(Ml) and gi(M„ 
(< (lMeV) 2 ). One can ignore k^M^), too. fi sub 



I) and gi(M 2 ) 
770MeV is what we have taken as 
the subtraction scale. In the case with L = T/2 = 2fm (where j3\ = 0.0836), one obtains 
£ = 1.3S and F = 1.2F, respectively. 

For the zero-mode integral, we use the partition function 



1 



x det 



2(^ 2 

-H b K v+ i{n b ) n v I v+ x(n v ) fil v+ i(fi) I v (ji) 

fllK u+2 (fJ, b ) ll 2 v I v+ 2{Hv) fJ 2 I u +2(fJ-) fJ,I v+ l(fJ,) 

\-^ b K u+3 (fi b ) nll v +z(n v ) /j?I u+3 (n) (i 2 I v+2 (ji) j 



(5.8) 



where [i = m u T,V = m^SV '. 

When the valence masses are degenerate, we use 



A„££ Q (/i„,/i) 



x£ Q (/i„,//) 



£ n v Y> 



Ji = 8m, 



£ 



£V 



(5.9) 
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We now present the explicit form of the correlators for the pseudoscalar and axial 
vector channels with m v = m v r, 



V c vv ,(t) = L 



?> 



2fi v £ 6F 2 M%V £ 



A c vv ,{t) 



_2^ S^(^ i /x) T/ti(t/r); (5n) 

where /ij = rrijSV '. We plot these correlators in Figs. || and |I] using &oo = 0.08331 for this 
case. 

5.2 iVi = 0, N h = 3 

As the second example, let us consider the case with m u = = 30 MeV, m s = 110 MeV 
while the valence quark masses are taken to be very light, m v = 0(1) MeV. In this case, 
all the sea quarks are in the p-regime and we therefore have iVj = and = 3. In this 
case, we have 

2{Ml d -Mlf _ 2{Ml d - Mlf Ml d Ml 

{Ml d + 2Mir l+ (Ml d + 2Mir M* d + 2Ml> ^ 

while M 2 = (Ml d + 2M 2 s )/3, where M 2 d = (m u + m d )S/F 2 , M 2 S = 2m s £/F 2 . Note that 
a double pole contribution now appears in r(t) as a partial quenching artifact since C ^ 0. 
The 1-loop corrections to the condensate and decay constant in this case are 



£ 1 

£ ~ 1_ F 2 



2Ml Ml ,„ », M 2 , M 2 , 

— i In + 2 gi (Ml) + — ^ In 
16?r M, u6 16tt 2 m, u6 



A / M 2 M 2 1 V „ ft , C f ln/^ 1 " , , 

" 3 I levr 2 111 ^ sub M 2 y I 3\/y 3 I 16vr 2 + P2 



-32^(/x su6 )(M 2 + M|;)j, 



F 1 
1 



F 2F 2 



2M 2 Ml , ~, Ml Ml w „ 9n 

— -f In + 2 51 M 2 + — ^ In - 16L^ Msu6 M 2 + M 2 - 

167r /4& 1Q7T Kub 



(5.13) 



where M 2 = (m u + m^)S/F 2 denotes the pion mass. Again we set = 770 MeV. In 
this case, the corrections are uncomfortably large: £ = 1.5£ and F = 1.2F. 
The zero-mode partition function for N[ = is given by 

Z iA + iN,=0)\m^) = det . (5.14) 
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When the valence masses are degenerate, we use 

/v '+ _ z v ' - ~~ „V ' +~ ^ „2> 

JT° = 2, .7 = 0, Ji=8m v 5^). (5.15) 
Here we present the correlators for the pseudoscalar and axial vector channels (m v = 

77V ), 



2fl v E 3F 2 M|y E i 

E 2 /2Cg^XV 



2F 2 3 E 



T 3 h 2 (t/T), (5.16) 



^(t) = — - ^^^)rfc!(t/r), (5.i7) 

where jli = rriiflV. We plot these correlators in Figs. [| and ^. 



6. Conclusions 

We have developed a new scheme of calculations for chiral perturbation theory with non- 
degenerate quark masses in finite volume. With our new counting rules separating iVj 
light quarks in the e-regime and the other Nh quarks in the p-regime, we have calculated 
the meson correlators in various channels: pseudoscalar, scalar, vector and axial vector. 
We have also calculated the disconnected contributions for the pseudoscalar and scalar 
channels. 

With the help of the replica method, we have also extended our study to the partially 
quenched case. Our results are shown to be consistent with all earlier work in the literature, 
both in the quenched and full QCD limit, with degenerate valence quarks. 

Our results can be compared to lattice QCD simulation with 2+1 flavors where the up 
and down quark masses are very light, but the volume is such that the theory is in the e- 
regime with respect to the corresponding pseudo Nambu-Goldstone bosons. The two-point 
functions are useful to determine the leading low-energy constants, the chiral condensate E, 
and the pion decay constant F in the chiral limit. As we have demonstrated, the formulae 
may also be used to extract the numerical values of higher-order low energy constants L\ 
and Lq. This work can be extended to other observables in the case where one valence 
quark is heavy, i.e. the chiral dynamics of kaons. With the new partially quenched chiral 
perturbation theory in this mixed regime one has an excellent analytical tool with which 
to explore future lattice simulations with nearly massless u and d quarks. It would be most 
interesting to investigate by analytic means also the region between the two regimes, where 
m n L ~ 1. 
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A. Space-integrals involving propagators 

In this appendix, we list several useful formulae for zero-momentum projection (or, equiv- 
alently, 3-dimensional space integrals) of functions expressed by A(x,M 2 ). 
A useful identity is 

^(2^) 2 + M 2 4Msinh(M^) L yV 

(A.l) 

which holds for an arbitrary regular function g(p). The zero-mode projection of A(x, M 2 ), 
for example, can be easily derived by setting g = 1; 

[ d 3 xA(x M 2 ) - 1 cosh(M(t - T/2)) _ _J_ 
J d XA(X ' M } ~ 2M sinh(MT/2) M 2 T' (A " 2j 

We are particularly interested in the massless limit: 

J d 3 xA(x, 0) = | (1 - I)' - L = Th^t/T). (A.3) 
It follows its second time derivative is given by 

fd"xd 2 A(x M 2 ) - M coshmt-T/2)) 
J d xd o^ M ) ~ 2 sinh(MT/2) ' 

J d 3 xd 2 A(x,0) = 1. (A.5) 

In this paper, we have also needed the following integral involving two A's 

J d 3 x J d 4 zd A(z - x,0)d A(z,0) = Th 1 (t/T), (A.6) 

and the more non-trivial integral 

J d 3 x (d A(x, M 2 )d A(x, M 2 ) - A(x, M 2 )8qA(x, M 2 )) 

- v ) + m2t2 ) + y y 2M sinh(MT/2) M 2 T J ' [A - () 
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where 



1 



, 4sinh 2 (v / |q| 2 + M 2 T/2) ' 
The chiral limit of Eq. ( |A.7| ) is given by 

d 3 x (d A(x, 0)<9 A(s, 0) - A(x, O)^ A(x, 0)) 



(Ai 



where 



A;, 



oo 



E 



|&oo+^i(vn 



-i 



i 

, 4sinh 2 (|q|T/2) + 12' 
becomes now a constant depending only on the shape of the box [ 



(A.9) 
(A.10) 



B. Summary of zero-mode group integrals 

Here we summarize the most essential zero-mode group integrals which are needed in the 



general partially quenched case, see also ref. [15| for additional details. 

The zero-mode contribution to the partition function with n bosons and m fermions 
is known as seen in Eq. fl4.ip . In this paper, we need the case with (n,m) = (1, N + 1) (N 
is the number of physical quarks): 



Z l,l+N{^b\^v,{^s}) 



nfi=i(^i - /4) nf 2 >s3(/ i s2 - k 



N 



x det 



si , 



LifJ-sz) ■ ■ ■ \ 



—HbK u+ i([Ib) ^vlu+l{^v) HslIu+l{Hsl) Ms2-fw+l(Ms2) 
filK v+2 (pb) ^llp+2{^v) ^1^+2(^1) ^2^+2(^2) 



...J 



(B.l) 



and (n,m) = (2,N + 2): 



^2,2+N 



j l^b2 \l^v\ 1 ^v2 1 

1 



0& 2 - /^i)(a4 - nfi=i(^i - a4)(a4 - nS> s3 (^2 - /4) 

/" Kv({l b i) Ku(nb2) Iu(fJ-vl) Iy{Pv2) ' ' \ 

—fl'blKv+lifl'bl) —V-b2K v +\(Hb2) Vvllv+li^vl) Mv2-Wl (^2) " 
Mfcl ^+2(^61) ^l 2 K v+2{^b2) nlxI v +2(Hv\) ^2^+2(^2) • • 



x det 



V 



/ 



(B.2) 



Here = mf,£V, = m„SV, where m&, m,, denote the masses of the valence bosons, 
the valence quarks respectively. Partially quenched observables can be computed by dif- 



ferentiating Eq. ( E^l ) or ( [B.2D with respect to suitable sources and subsequently taking the 
limit fif, — > fj, v . 
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As building blocks, we use two quantities denned in Eqs. (|4.3| ) and ([4^ 

^ Q( ^' {M) and I^ Q (/W,2,{M). (B.3) 
Note that in the degenerate limit pL v \ = fj, V 2 = fJ* v , 

(B.4) 

The needed formulae for one valence index are 

.PQ 



2 ((A™ + Ul v ))u = , (B.5) 

l /m _ L ^rt^2\ d M ^ Q (/^,{M) AS^ Q (^,{^ S }) 

= (B.7) 

1 HTT TT\ \ 2 \ — S ^ Q (^' Off}) i ^ 2 /"r> o\ 



(A«C4,>1/ = |((A« + Ulf)u - \{{U VV - Ulf)u 



a M „S^ Q (/i„,{/i 8 » _ As£ Q (/i„,{/x fl » ££ Q (/x„,{/A a }) _ i/* 

For two valence indices, 

~((Z7 OT + Ul){XJ v , v , + Ul, v ,))u = I>P Q (^,AV, W)i (B.10) 

~<(17 W - t&KEW " = 7^T> (B-ll) 

(U vv U v , v ,)u + {Ul v Ul, v ,)u = 2^0^, /v,M) + • (B.12) 

Similarly, 

\{{U VV , ± 17+ J2 )p = l ((f/ ^ ± ^2^ 



±1 



x£ Q (/i„,{M) E^ Q (/v,{M) \ 
^ Mo' = 1 



\(Ulr + (Ul, v f)u = 0, (B.13) 
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as well as 



\{{u m , ±ul, v )(jj v , v ±ul,)) u = 

1 / S^ Q (/x„,{/z 8 » ^%v^})\ 

/V ^ ^ , (B.14) 



..2 ,,2 W v f*« y; 

Vv V v ' \ ^ ^ / 

(tW^ + U^U^u = 0. (B.15 



were also derived in ref. [15|. 



C. Some Ward-Takahashi Identities at fixed topology 

For the computation of vector and axial vector correlation functions we needed a set of 
zero-mode expectation values involving three zero-mode fields U. These can be reduced 
to known integrals by means of exact identities on the group manifold of U(N). Such 
relations correspond to Schwinger-Dyson equations on the group manifold and encode, in 
physics terms, Ward-Takahashi Identities (WTI) of spontaneous chiral symmetry breaking 
in a sector of fixed topological charge v. The derivation below follows the method described 
in detail in Appendix B of ref. pJ] . 

Let t a denote generators of U(N) in a chosen representation, here the fundamental. In 
addition, let e a be infinitesimal parameters. We introduce left-handed differentiation V a 
on the group by means of 

F(e ieata U) = F{U) + e a V a F(U) + . . . (C.l) 

The derivatives V a give rise to a standard Leibniz rule, and left-invariance of the Haar 
measure on U (N) ensures that 

J dU V a F(U) = . (C.2) 

Choosing different functions F{U) this simple identity generates an infinity of exact rela- 
tions on the coset of symmetry breaking for the zero-mode fields. For the present purposes 
we can choose, e.g., 

F{U) = Tr[M 1 [/]Tr[[/ t M 2 ]P(?7), (C.3) 

where M\ and M 2 are arbitrary N xN matrices, and the Boltzmann weight P(U) is defined 
in the obvious way: 



P(U) = (detC/yexp 



—tr(MU + U ] M) 



(C.4) 



Different choices of the matrices M\ and M 2 lead to identities that are useful in con- 
nection with the vector and axial vector correlators. For example, (Mi)y- = Si v it",j and 
(M 2 )ij = Si v 5 v j (and the similar choice with indices v and v' swapped) gives, after use of 
the U(N) completeness relation (with a sum over a) (t a )ij(t a )ki = ^SySjk, 

(Ul(UMU) vlvl )u = (m v > v >UI - ^(N + v)UlU v , v ^ , (C.5) 
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where we used = M , and the hermitian conjug ate relation (Note (det U) u = (det U^)~ v ), 

(U vv (U^MU%, v ,)u = (m vV U w - ^{N - v)U vv ul, v )j . (C.6) 

Another choice, {Mx)%j = ^t^,. and (M 2 )ij = 5 iv 8 v ij gives 

(UI, V (UMU) V/V )u = ~(N + v) (ul v U v , v ) v = 0. (C.7) 
See the appendix [B] for the last equality to zero. The hermitian conjugate is also vanishes, 

(U v , v (U^MU% v }u = 0. (C.8) 
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v=0, m=2 MeV, m=l 10 MeV 




1 1 1 1 1 — 1 

0.2 0.4 0.6 0.8 1 

t/T 

Figure 3: The pseudoscalar correlators with m v = 1-3 MeV, m u = = 2 MeV and m s = 110 
MeV in a sector of trivial topology, v — (top) and in sectors of v = 0-2 for fixed m v = 1 MeV 
(bottom). We set L = T/2 = 2 fm and the correlators are normalized by the Sommer scale ro = 0.49 
fm §§. 
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v=0, m u =2 MeV, m=l 10 MeV 
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Figure 4: The axial correlators for same parameter values as in Fig. ||. 



v=0, m M =30 MeV, m=\\0 MeV 
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Figure 5: The pseudoscalar correlators with m v = 1-3 MeV, m u = nrid = 30 MeV and m s = 110 
MeV in a sector of trivial topology, v = (top) and in sectors of v = 0-2 for fixed m v — 3 MeV 
(bottom). 



- 31 - 



v=0, m M =30 MeV, m=UO MeV 
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Figure 6: The axial correlators for the same parameter values as in Fig. pj. 



